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Systems of strongly correlated fermions on certain geometrically frustrated lattices at particular 
filling factors support excitations with fractional charges zbe/2. We calculate quantum mechanical 
ground states, low-lying excitations and spectral functions of finite lattices by means of numerical 
diagonalization. The ground state of the most thoroughfully studied case, the criss-crossed checker- 
board lattice, is degenerate and shows long-range order. Static fractional charges are confined by 
a weak linear force, most probably leading to bound states of large spatial extent. Consequently, 
the quasi-particle weight is reduced, which reflects the internal dynamics of the fractionally charged 
excitations. By using an additional parameter, we fine-tune the system to a special point at which 
fractional charges are manifestly deconfined — the so-called Rokhsar-Kivelson point. For a deeper 
understanding of the low-energy physics of these models and for numerical advantages, several 
conserved quantum numbers are identified. 



I. INTRODUCTION 

Quantization of charge is a very basic feature in the description of the physical world. Therefore, 
fractionally charged excitations came as a surprise to physicists. Already back in the year 1979, Su, 
SchriefTer, and Heeger 1 showed that a model describing the one-dimensional (ID) chain molecule 
trans-polyacetylene (CH) n supports excitations with spin-charge separation. This is not yet charge 
fractionalization, but when the model is considered at different electron densities (corresponding 
to extremely high doping), it turns out that at certain filling factors elementary excitations with 
fractional charge exist [in the simplest case q = d=e/3 and q = d=2e/3]. A few years later, Laughlin 2 
interpreted the much celebrated fractional quantum Hall effect (FQHE) in terms of fractionally 
charged (quasi-) particles (fcp) and fractionally charged (quasi-) holes (fch). Thereby, he firmly 
established the idea of fractional charges in our understanding of solid-state physics. Direct electron- 
electron interactions are crucial in the FQHE: Correlations become strong in an applied magnetic 
field because the kinetic energy of the electrons is quenched. Of course, no one will ever extract a 
fraction of an electron from an FQHE sample. But thinking of an extra electron or an extra hole near 
fill factor v = 1/3 as decaying into three separate entities of charge q = d=e/3 each proved enormously 
helpful for a qualitative and quantitative understanding of the FQHE — and was rewarded by the 
Nobel Prize in Physics 1998. 

The question was left open whether or not fractionally charged excitations exist in 2D or 3D 
systems without a magnetic field. In 2002, it was suggested by one of us^ that excitations with 
charge ±e/2 do exist in certain 2D and 3D lattices, e.g., the pyrochlore lattice, which is a prototype 
3D structure with geometrical frustration. The original work was motivated by the transition metal 
compound LiV204, which surprisingly shows heavy-fermion behavior with, e.g., a large 7 coefficient 
in the low-temperature specific heat C = 7T. 4 However, we will not address the issue whether the 
models discussed here apply to any particular material or artifical systems such as optical lattices, see 
e.g. Ref£ and citations therein. Instead, we try to contribute to the very general question whether 
or not fractional charges can exist at all in truly 2D or 3D systems in the absence of magnetic fields. 

We would like to argue that general prerequisites for fractional charges are — as in the FQHE case — 
strong short-range correlations and certain band fillings. Furthermore, the short-range correlations 
should be somehow incompatible with the lattice structure in order to prevent the development 
of long-range order. Following the general usage, we will simply call the lattice "frustrated" even 
though calling the interaction "frustrated by the lattice" would be the more accurate terminology. 

In this contribution, we study the charge degrees of freedom on such lattices systematically and 
consider a class of models of strongly correlated spinless fermions. Most of our numerical calculations 
were done for the 2D checkerboard lattice, which is easier to deal with than the even more interesting 
3D cases. For future work, one can hope to learn from comparison with spin systems where numerous 
studies exist for the pyrochlore lattice and its 2D relatives, the crisscrossed checkerboard lattice and 



Figure 1: (a) Checkerboard lattice as 2D projection of the 3D pyrochlore lattice.— (b,c) Two examples of 
allowed configurations on a checkerboard lattice at half filling. Occupied sites are connected by thick lines 
as guides to the eye. 



the kagome lattice. 
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II. FRACTIONAL CHARGES ON FRUSTRATED LATTICES 



In order to illustrate the concept of fractional charges on frustrated lattices, we focus here on the 
(crisscrossed) 2D checkerboard lattice. Figure QJa) illustrates that the checkerboard lattice can be 
thought of as a projection of the pyrochlore lattice onto a plane. In the following, we adopt the 
ideas of Ref. 3 and consider a model Hamiltonian of spinless fermions 



on a checkerboard lattice. The operators c^(q) annihilate (create) fermions on sites i. The density 
operators are ni = c\c i . We assume half filling, ni = | for a system with N sites. Our 
main interest is the regime where the nearest-neighbor hopping t is much smaller than the nearest- 
neighbor repulsion V, i.e., \t\/V <C 1. Henceforth, we assume t > 0. In analogy to the tetrahedra in 
a pyrochlore lattice, all bonds in a crossed square are equivalent. 

a. Classical correlations and ground-state degeneracy. For a moment, let us set the hopping- 
matrix element t to zero. The ground-state manifold is then macros copically degenerate: Every 
configuration that satisfies the so-called tetrahedron rule of having exactly two particles on each 
tetrahedron (crisscrossed square) is a (classical) ground stated and will henceforth be referred to 
as an "allowed configuration." Examples are shown in Fig. [2(b) and (c). Note that our coordinate 
system is rotated by 45° relative to that of, e.g., Ref. 13 . The resulting difference in boundary 
conditions can lead to noticeable numerical differences in particular for small cluster sizes. 

One can visualize the origin of the macroscopic degeneracy as follows: Take the set of the six 
allowed crisscrossed squares and construct row by row a larger allowed configuration. Whenever a 
crisscrossed square is added, we can choose between one, two, or three different possibilities, depend- 
ing on the neighboring crisscrossed squares. Since there often is a choice to make, an exponential 
number of different allowed configurations can be constructed. Thus, the system has a finite entropy 
at T = 0. The exact value of the ground- state degeneracy can be obtained from a mapping to the 
so-called six- vertex model. 14 Its solution is highly non-trivial due to the existence of long-range 
correlations, which seem to be a generic feature seen in many frustrated lattice models. 15 

All classical (t = 0) ground states have the important property of being incompressible in the 
sense that no fermion can hop to another empty site without creating defects and thereby increasing 
the repulsion energy. In other words, we have to violate the tetrahedron rule in intermediate states 
if we want to connect via hopping processes one allowed configuration with another. 

b. Fractional charges. Placing one additional particle with charge e onto an empty site of an al- 
lowed configuration leads to a violation of the tetrahedron rule on two adjacent crisscrossed squares, 
see Fig. Efa). The energy is increased by AV since the added particle has four nearest neighbors 
(charge gap). There is no way to remove the violations of the tetrahedron rule by moving electrons. 
However, fermions on a crisscrossed square with three particles can hop to another neighboring criss- 
crossed square without creating additional violations of the tetrahedron rule, i.e., without increase 
of the repulsion energy [see Fig. [2](b,c)]. By these hopping processes, two local defects (violations 
of the tetrahedron rule) can separate and the added fermion with charge e breaks into two pieces. 




(i) 
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Figure 2: (a) Adding one fermion to the half-filled checkerboard lattice leads to two defects (marked by 
black triangles) on adjacent crisscrossed squares, (b), (c) Two defects with charge e/2 can separate without 
creating additional defects. They are connected by a string consisting of an odd number of fermions. 




Figure 3: (a)-(c) Hopping of a fermion to a neighboring site in an allowed configuration generates a fluctu- 
ation: (a) A fractional charged particle (fcp) and fractional charged hole (fch) are generated, (b)-(c) The 
two defect (marked by triangles) with charge zbe/2 can separate without creating additional defects and are 
connected by a string consisting of an even number of fermions. (d) Example of an allowed configuration on 
a checkerboard lattice at half filling with possible low-order hopping processes. 

They carry a fractional charge of e/2 each. In the quantum mechanical case (t ^ 0), the separation 
leads to a lowering of the kinetic energy of order \t\. 

Energy and momentum must be conserved by the decay processes le — > 2fcp / s. If we associate 
momentum k and energy E(k) with the inserted fermion, they must now be shared between the 
resulting two fcp's into which it has decayed 

E(k)=W + e(k 1 )+e(k 2 ). (2) 

Here k = ki+k2 and e (k) is the energy dispersion of a fcp. Even though e (k) is at present completely 
unknown, Eq. ([2j) allows to predict that for deconfined fcp's the electronic spectral function should 
not contain a Fermi-liquid peak, but should show a broad continuum instead. 

Figure [2] demonstrates that the two defects can be thought of as always being connected by a 
string of occupied sites consisting of an odd number of sites. The fractional charges can thus be 
alternatively interpreted as the ends of a string-like excitation. In this picture, the connection is 
not static, as two pairs of fcp's can exchanges partners when the connecting strings come close. 

c. Quantum Fluctuations. If we relax the constraint of having two fermions on each crisscrossed 
square and consider a small but finite ratio t/V, quantum fluctuations come into play. The quan- 
tum fluctuations lead also to fractional charges, but do not change the net charge of the systems. 
Starting from an allowed configuration, the hopping of a fermion to a neighboring site increases the 
energy by V. One crisscrossed square contains three fermions while the other has only one fermion, 
see Fig. E^a)-(c). These so-called vacuum fluctuations (virtual fcp-fch pairs) lead to two mobile 
fractional charges with opposite charges -fe/2 and —e/2, which are connected by a string of an even 
number of fermions. The energy associated with a free fcp and a free fch is A£? V ac = V+e (k)+e (— k), 
where e (— k) denotes the kinetic energy of a fch. Such virtual processes connect different allowed 
configurations, lower the total energy, and reduce the macroscopic degeneracy, as we will see more 
explicitely in the next section. 

d. Fractional charges in 3D. All arguments mentioned above for the existence of fcp's on a 2D 
checkerboard lattice can be directly transferred to a 3D pyrochlore lattice at half filling, see Fig. 2J 
Fcp's correspond to tetrahedra with three fermions and fch's to tetrahedra with only one fermion. 
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(a) (b) 

Figure 4: Adding one fermion to the half-filled three-dimensional pyrochlore lattice leads to two defects on 
adjacent tetrahedra. The two defects with charge e/2 can separate without creating additional defects and 
are connected by a string consisting of an odd number of fermions. 

III. EFFECTIVE HAMILTONIANS 

e. Ring exchange processes in the undoped case. The numerical and analytical work is greatly 
simplified, if a simpler Hamiltonian H can be derived that shows in the limit |t| <C V the same low- 
energy physics as the model Hamiltonian (pQ). A down-folding procedure can be used to define TL\ Per 
definition, it acts only on the subspace of allowed configurations and includes all virtual processes up 
to the lowest non-trivial order, which is t 3 /V 2 . The relevant processes are shown in see Fig.[3fd) and 
can be classified as self-energy contributions and ring-exchange processes i^ e ff, W = H% + H e ^ . 
The former comprises the terms which are diagonal in the real space basis. The latter includes those 
which connect different allowed configurations. H% contains fermion hops to an empty neighboring 
site and back again as well as hops around an adjacent triangle. These contribute only a constant, 
configuration-independent energy shift. It does not lift the macroscopic degeneracy and will hence 
be ignored furtheron. The total amplitude of ring-exchange processes around empty squares is 
proportional to t 2 /V. It vanishes for spinless fermions because the amplitudes for clockwise and 
counter-clockwise ring-exchange cancel each other due to fermionic anti-commutation relations. 
Thus, the macroscopic ground-state degeneracy is first lifted by ring exchanges ~ t 3 /V 2 around 
hexagons, and the looked-for effective Hamiltonian reads 

ffeff = -5E {o , }(l < X > >«^ > |-| < X>>«^ > l +H - C -)' ( 3 ) 

with the effective hopping-matrix element g = 12 t 3 /V 2 > and the sum taken over all vertical and 
horizontal oriented hexagons. The pictographic operators represent the hopping around hexagons 
which have either an empty or an occupied central site. The signs of the matrix elements depend on 
the representation and the sequence in which the fermions are ordered. When the sites are enumer- 
ated along diagonal rows, an exchange process commutes an odd number of fermionic operators if 
the site in the center of the hexagon is empty and an even number if the center is occupied. In Ref. 13 , 
it has been shown that the effective Hamiltonian (|3]) gives a good approximation of the low-energy 
excitations of the full Hamiltonian (pQ) in the limit considered. 

/. Propagation of defects and extra charges. We have derived H e ^ as effective Hamiltonian for 
local rearrangement processes as resulting from the (virtual/intermediate) generation of a defect pair 
which subsequently recombines, leaving behind a different allowed generation. For many questions 
of physical interest, it is necessary to consider the propagation of defects over large distances. Anal- 
ogously to conventional semiconductors, this refers to long-lived thermally generated defects (e-h 
pairs in the semiconductor analogy) as well as to the consequences of slight doping, i.e., addition of 
two fcp's by addition of one extra fermion. The natural generalization of the effective Hamiltonian 
(|3j) to these cases includes besides the lowest order ring-exchange processes H e fi a projected hopping 
term that moves the defects 

#tg = H eE -tY,p + H - c -) R ( 4 ) 

(id) 

The projector P ensures that H tg acts only on the subspace of configurations with the smallest pos- 
sible number of violations of the tetrahedron rule which is compatible with the number of particles, 



5 




Figure 5: (a) Height representation for examples of allowed configurations of a V32 x checkerboard 
lattice with periodic boundary conditions at half filling. The height field h (numbers in the non-crossed 
squares) is uniquely defined for a given configuration up to an additive constant M. The field f = V/i is 
indicated by small arrows on the lattice sites. Details of the mapping can be found in the text. The effect 
on the height fields of a ring-exchange process around a hexagon is shown explicitely. (b) The effective 
Hamiltonian conserves the number of fermions on each of the four sublattices of the checkerboard lattice 
which are labeled by "blue squares" , "red circles" , "yellow stars" and "green triangles." (c) Ring-exchange 
processes change the number of fermions on sites marked by dark "purple squares" by two. 16 

e.g., two in the case of one added fermion. We refer to H tg as the t-g model and consider the param- 
eters t and g as independent (i.e., not restricted to the regime g <C t as enforced by t < 7). This 
allows to study the effect of ring exchange onto the dynamics of fractionally charged excitations. In 
particular, the question of confinement can be studied even on rather small clusters by increasing 
the ring exchange strength g relative to t. 

IV. HEIGHT REPRESENTATION, CONSERVED QUANTITIES, AND GAUGE 

SYMMETRIES 

Conserved quantities of a Hamiltonian allow for a reduction of the numerical effort by exploiting 
the resulting block-diagonal form of the matrix representation. The real space configurations span- 
ning the subspaces corresponding to these blocks will be referred to as "subensembles." Eigenstates 
can conveniently be classified by the eigenvalues of the conserved quantities as quantum numbers 
for our model. We will now identify several such quantum numbers in order to characterize different 
subensembles. A topological quantity, which is conserved by all local processes, i.e., ring-exchange 
processes, is the average tilt of a scalar height field which will be introduced in the next paragraph. 
Another useful set of quantum numbers are the number Ny, Ng, Nn) of particles on the four 
sublattices referred to as blue, yellow, green and red as shown in Fig. 03(b). They are conserved by 

Any allowed configurations of the half-filled checkerboard can uniquely be represented by a vector 
field f for which the discretized lattice version of the (discrete) curl vanishes, i.e., a pure (discrete) 
gradient of a scalar field (height field) f = Vft, see Fig. [5] The height field h, is derived from 
the local constraint expressed by the tetrahedron rule as follows: A clockwise or counter-clockwise 
orientation is assigned alternatingly to the crisscrossed squares. Arrows of unit length are placed 
on the lattice sites. The arrows point along (against) the orientation of the adjacent crisscrossed 
squares if the site is occupied (empty) . It is easily checked that allowed configurations are those for 
which the discretized line integral of f around every closed loop vanishes. Thus, f = Vft defines a 
height field h up to an arbitrary constant M. 

The height at the upper (right) and at the lower (left) boundary of a finite lattice of N x x N y 
squares with periodic boundary conditions can differ only by an integer —N y ^ < K y ( x ) < N y ( x )^ 
which is the same for all columns (rows). This defines topological quantum numbers (n Xl n y ). They 
remain unchanged by all local processes that transform one allowed configuration into another, i.e., by 
ring-exchange processes along contractible loops. In particular, H e R merely lowers or raises the local 
height of two adjacent plain squares by ±2, as illustrated in Fig. [5](a). We refer to (k x /N x , K y /N y ) 
as global slope. 

The lattice symmetry is broken for states with (k x , n y ) ^ (0,0): 13 E.g., a finite positive value of 
k>x/N x implies a charge modulation along diagonal stripes. Similarly, a charge density modulation 
is present if the condition N& = Ny = Ng = Nn is violated. 

For the identification of exactly solvable points in parameter space and for future quantum Monte 
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Figure 6: The three different actions of the effective Hamiltonian on the topology of loop configurations are 
shown in panel (a). Panels (b) and (c) show representations of two configurations by fully-packed directed 
loops.— 

Carlo simulations, it would be very advantageous if all matrix elements had the same sign. As first 
steps in this direction, we describe a gauge transformation that changes the global sign of g in the 
effective Hamiltonian ([3]) and then show conditions under which it is possible to remove minus signs 
for the half-filled checkerboard latticed 

Consider a "purple" sublattice V that contains the sites around every second plaquette as shown 
in Fig. Etc). Define cr-p as the number of fermions on purple sites for a given configuration: <j<p = 
Y^iev Ui ' One observes that ring-exchange processes change a-p by two. Thus, if all configurations 
are multiplied by a factor of i av , the sign of all ring-exchange matrix-elements are changed. The 
invariance with respect to this gauge transformation proves a global g <-> (—g) symmetry. 

However, a fermionic sign problem remains: Ring-exchange processes around empty and occupied 
hexagons carry opposite signs in Eq. ([3|). We argue that it can be avoided in certain (but not all) 
cases iii In order to do so, we represent the ground-state manifold by ensembles of fully-packed 
loops as exemplified in Fig. dJ We notice that a ring exchange around a hexagon with an occupied 
center site does not change the loop topology, whereas a ring exchange around an empty hexagon 
always does cause changes in one of the three topological different ways shown in Fig. [61(a). 

Let us consider allowed configurations with "fixed" boundary conditions with an even number of 
fermions on the four boundaries. They are represented by closed loops in the interior and loops 
terminating at a boundary. We orient the closed loops as follows: (i) Color the areas separated by 
the loops alternatively white and grey, with white being the outmost color; (ii) orient all loops so that 
the white regions are always to the right, see Fig. [BJb). If open loops are present [Fig. 15(c)], these 
are closed arbitrarily but intersection- free outside the sample and colored as described above. Let 
us assume without loss of generality that the color at infinity is white. We now notice by inspection 
of Fig. ([6]) that the relative signs resulting from the exchange processes around empty hexagons are 
consistent with multiplying each loop configuration by i r (— i) 1 , where r and I are the total number of 
the clockwise and counter-clockwise winding loops, respectively. Hence, by simultaneously changing 
the sign of the exchange-processes around empty hexagons and transforming the loop states 

\C) - i l ^(-i) r ^\C), (5) 

we cure the sign problem, thus making the system effectively bosonic. 

This construction need not work for periodic boundary conditions: Firstly, only even-winding 
subensembles (sectors) on a torus allow for such a two-color coloring. Secondly, even then it might 
be possible to dynamically reverse the coloring while returning to the same loop configuration. 
However, the exact diagonalization results presented below (see Fig. [8j) suggest that for periodic 
boundary conditions on even tori (preserving the bipartiteness of the lattice), the lowest-energy 
states belong to a sector where such a transformation works. We remark that the presented non-local 
loop-orienting construction is restricted to the effective Hamiltonian ([3]), i.e., to the ring-exchange 
processes of length six. 



Figure 7: (a)-(b) Fragments of possible "frozen" ground states (no flippable hexagons), (c) One of several 
configurations of the half-filled checkerboard lattice that maximize the number of flippable hexagons. The 
unit cell contains \^20 x \^20 sites. 



V. GROUND STATES AND LOWEST EXCITATIONS IN THE UNDOPED CASE 

In order to discuss the possible confinement of fcp's, we have to investigate the nature of the 
quantum-mechanical ground state of the undoped system. We do so in the approximation of an 
effective ring-exchange Hamiltonian ([3|) acting only on the subset of allowed configurations. Following 
Rokhsar and Kivelson^ we add to H e ^ an extra term that counts the number of flippable hexagons. 
The extended Hamiltonian reads 

H gll = ^e ff + ^E {o ,0}(l<X>)(<^>l + l<X>)«X>l). ( 6 ) 

where the pictographic operators with grey-colored dots in the center are summed over all flippable 
hexagons, independent of the occupancy of the site in the center. Next, we discuss some limiting 
cases of the Hamiltonian (|6]): 

(i) \i — > +00: All configurations which contain no flippable hexagons (frozen configurations) become 
ground states with E = 0. Some are shown as Fig. [71 (a,b). 

(ii) \i — > — 00: Ground states are configurations with maximal number of flippable hexagons iVfl. 
Using a simple Metropolis-like Monte Carlo algorithm for lattices with up to 1000 sites, we always 
find that configurations of the type shown in Fig. [71(c) or slight variations thereof maximize Nr. Such 
configurations will be referred to as "squiggle" configurations. 19 Our numerical results suggest that 
in the thermodynamic limit the degenerate ground-state configurations all lie in the (k x , k v ) = (0, 0) 
subspace. 

(hi) \i = g > 0: This is the exactly solvable Rokhsar-Kivelson (RK) point. 18 Following the original 
RK construction, we rewrite the Hamiltonian ([6]) for \i = g in a way which explictly shows that some 
liquid-like ground states have energy E = and, thus, become degenerate with the frozen states, 
i.e., the solutions. We assume that we can use the gauge transformation ([5]) to change the 

sign in the second term and to rewrite the Hamiltonian as 

H g =, = 9 E [(l<X»-KS»)x((<S>|-«X>l)]- (7) 

Since this is a sum over projectors, all eigenvalues are non-negative. Furthermore, after re- 
gauging, all off-diagonal elements are non-positive and for each subensemble t an exact ground- 
state wavefunction is given by the equally weighted superposition of its configurations |C^), i.e., 
\^rk) ~ Si I^P)- These coherent superpositions are the analogs of the Resonating Valence Bond 
(RVB) state, originally discussed by L. Pauling 20 and P. W. Anderson* 2 ^ Their energy is easily 

computed: Each flippable hexagon contributes — g and -h/i, thus for g = \i\ {^q^k^z^q^rk) ~ 0- 
Note that for fermionic systems, we find a well defined RK point only if a gauge transformation 
exists such that all off-diagonal matrix elements are non-positive. Otherwise the energy is most 
likely larger than zero and the subensembles consequently do not form a ground state. 

Next, we explore by means of numerically exact diagonalization the eigenstates of small clusters 
general /i values. For the actual calculations, we first generate all configurations that fulfill the tetra- 
hedron rule, then group them according to quantum numbers and generate a sparse block-diagonal 
matrix representation of the Hamiltonian (|6|). For a 72-site checkerboard cluster with periodic 
boundary conditions, the 16 448 400 dimensional low-energy Hilbert space of allowed configurations 
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Figure 8: (a) Energies of the ground state and lowest excited states in each subensemble of a 72-site half- 
filled checkerboard cluster for different values of /x of the g-\i Hamiltonian. Level crossing of ground states 
occurs only at ji = g. The insets indicate different phases: Maximal flippable plus fluctuations for \i < g, 
a critical point \i — g where the ground state is an equally weighted superposition of all configurations, 
and frozen configurations as ground states for fi > g. (b) Left side: Ground-state energy and energies of 
the lowest excited states of the effective Hamiltonian if e ff in subspaces with different global slopes (K x ,K y ) 
= (k x ,0) for 72-site cluster. Right side: Same system, but assuming same signs for all matrix elements 
("bosonic calculation"). 17 



can be decomposed into a few hundred subspaces, where the largest one has 1 211 016 dimensions. 
The low-energy states of the sparse block-diagonal Hamiltonian matrix and physical properties such 
as charge density distribution and density-density correlations can easily be obtained on a 64-bit 
workstation. Figure [5(a) shows energies of the ground-state and the lowest excited states of all 
subensembles. 

(iv) For fi < g, two ground states are found in subensembles with (k x ,k v ) = (0,0) and 
{N& , Ny , Ng , Nfi) = (6,6,12,12) and (12,12,6,6), respectively. These are superpositions of con- 
figurations with the maximal number of flippable hexagons. At the physical point \i = 0, the 72-site 
system is in a crystalline and confining phase. In the thermodynamic limit, we expect to recover 
the 10-fold degeneracy of the squiggle phase instead of the two-fold degeneracy. For small values 
H > 0, the average weight of configurations in the ground state with the maximum number iVfl is 
large (not shown), but decreases with increasing \i until at the RK point, i.e., /i=g, the ground state 
is formed by an equally weighted superposition of all configurations of a certain subensemble. 

(v) For fi > we find essentially the same properties as described above in the limit \i — > oo. 

In summary, we find a confining phase and a phase in which the ground states are given by static 
isolated configurations. The two phases are separated by a point with deconfined excitations, i.e., 
the Rokhsar-Kivelson (RK) point ^ The main finding is that the original effective Hamiltonian is in 
a confining phase with a long-range ordered ground state (squiggle phase). This phase maximizes 
the gain in kinetic energy and is stabilized by quantum fluctuations (order from disorder). The 
results of the exact diagonalization on small samples indicate that the symmetry remains broken 
all the way along the /i-axis up to the RK point. This observation is also strongly disfavoring a 
deconfining phase to the left of the RK point. Hence, the fermionic RK point is likely to be an 
isolated quantum critical point just as it is for the bosonic model 

We come back to the fermionic sign problem. Figure [8(b) compares energies of a system in 
which the Fermi sign is taken into account to those from calculations which exclude the Fermi sign. 
The ground-state energy, the first excited states in the (/%, n y ) = (0,0) sector, and the weights of 
the different configurations in the corresponding eigenstates are the same in both case. In some 
subensembles with (k x ,k v ) ^ (0,0) the energies of the ground states including the Fermi sign are 
higher than the ones for bosons (no Fermi sign). 

At the physical point \i = 0, the 72-site ground state is two-fold degenerate with quantum numbers 
(iVg, Ny, Ng, Nfc) = (6, 6, 12, 12) or (12, 12, 6, 6). The resulting charge order with alternating stripes 
of average occupation 1/3 and 2/3 is shown in Fig. [9(a). The density-density correlation function 
in the quantum-mechanical ground states l^o^) [Fig. [9(b)] 

is best understood as reflecting the algebraic correlations present in the average over all classically 



Figure 9: Half filling: (a) Charge density distribution for one of the two ground sates, (b) Corresponding 
density-density correlation function Cf^. The site io with average density 2/3 shows up as the largest dot 
in the panels (b)-(c). The radius of the dots is proportional to the absolute value. A red or blue color 
represents a positive or negative value, respectively, (c) Classical density-density correlation function. 
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Figure 10: (a), (b), and (e): Local loss of kinetic energy due to the separation of two fractionally (static) 
charged defects (fcp's marked by light red squares or fch's marked by dark blue squares). The radii of 
the circles are proportional to the local energy loss, (c), (d), and (f): Red (blue) circles show an increase 
(decrease) of the local density (vacuum polarization due to the two fcp's or the fcp-fch pair). 

degenerate configurations, see Fig. ^c) and Ref. 15 . Specific quantum-mechanical features resulting 
from ring hopping become visible in the difference of the actual correlations and have been discussed 
in Ref,ia. 



VI. STATIC AND DYNAMICAL PROPERTIES OF THE DOPED SYSTEM 



Next, we turn to the important issue of confinement /deconfinement of defects generated either 
particle-hole excitations or by weak doping. Two different kinds of calculations can be done: In 
the first case — referred to as static charges — the defects are fixed to some lattice positions and 
ground-state calculations involving H e R are performed in the restricted subspace spanned by the 
configurations with exactly those defects at given positions. The total (free) energy is recorded 
as function of the defect positions. Its spatial derivative is interpreted as attractive or repulsive 
force. In the second case — referred to as dynamic charges — the defects propagate according to the 
Hamiltonian (jl]). The existence or non-existence of bound states is interpreted as confinement or 
deconfinement. 

g. Static charges. If an additional particle is added/doped to a half-filled checkerboard lattice, 
at least two crisscrossed squares are occupied by three particles. It is easy to see that, e.g., for the 
t-fi Hamiltonian in the limit of large positive \i (frozen configurations) the ground-state energy of 
two static fractional charges is independent of the distance between. Thus static charges are not 
attracted by a force, which suggests that dynamic fcp's are deconfined at zero temperature in that 
limit. 
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Also, it is obvious that at the RK point the total energy of a system with two static charges is 
independent of the distance between them. Thus dynamic fcp's are expected to be deconfined. 

We will not discuss systems with a general values of n further, but focus on the physical case n = 
and discuss the changes of the kinetic energy density in the presence of two static charges 
The energy change can be decomposed into local contributions from all ring-exchange processes 
involving a given site i. An increase in kinetic energy in the region between the two fractional charges, 
i.e., along the connecting string, is found and illustrated in Fig. [TOla.b). The total energy increase 
is approximately proportional to the length of the generated string and implies at large distances a 
constant confining force. The changes of the local energy density goes along with density changes, 
as illustrated in Fig. [TOlc.d). Similar results are found when a particle is removed from the ground 
state or when a particle-hole excitation is generated out of the ground state, see Fig. [T0lfc,d,f). 
The calculations have been performed in a reduced Hilbert space in which only two static defects are 
present. In a calculation within the full Hilbert space, or a subspace allowing at least an additional 
fcp-fch pair, the energy would not increase linearly to infinity, but the connecting string is expected 
to break by creating additional pairs of defects when this energetically favorable. For the relevant 
parameters, e.g., g = O.Olt and V = lOt, this occurs when two fcp's have separated over 1000 
lattice sites. This effect is well known for the case of confined quarks where pair production (quark 
anti-quark pairs) occurs before the quarks have been separated to an observable distance. 

The attractive constant force acting between two fractional charges in the confining phase results 
from a reduction of vacuum fluctuations and a polarization of the vacuum in the vicinity of the 
connecting strings. These findings suggest that a number of features known from QCD are also 
expected to occur in a modified form in solid-state physics. Conversely, one would hope that by 
studying frustrated lattices or dimer models one might be able to obtain better insight into certain 
aspects of QCD. 

h. Dynamic charges. We turn now to the dynamical properties of fcp's, in particular spectral 
functions and optical conductivity for the half-filled checkerboard lattice. Numerical studies of 
finite clusters are again the method of choice^ because conventional approximation schemes such 
as mean-field theories or Green's function decoupling schemes are unable to describe the strong local 
correlations expressed by the tetrahedron rule. Diagonalizations with up to 50 sites were performed 
within the minimal Hilbert space spanned by the configurations with the smallest possible number 
of violations of the tetrahedron rule [half-filled system: all allowed configurations plus those with 
one fcp-fch pair; doped case: 2 fcp's or 2 fch's]. 

The spectral function A(k, co) = A~(k, co) + A + (k, co) of an interacting many-particle system is 
the sum of the probability amplitudes for adding (+) to the iV-particle ground-state system l^o^) 
or removing (-) from it a particle with momentum k and energy co (h = 1) 

A+(k, W ) = lim -llm(^|c k 1 4 (9) 

77^0+ 7T CO + IT] + h/Q — 12 

A"(k,o;) = lim ■ -jWffl 4 1 c k |<). (10) 

77^0+ 7T CO + IT] — £L/Q + £1 

We use operators c\. — Y^j e zr?k a| in the extended Brillouin zone. 

The spectral functions yield direct insight into the dynamics of a many-body system, as seen, e.g., 
in angular-resolved photoemission spectroscopy (ARPES). Expectation values of the form G(z) = 
(ipo\A (z — H)~ x At |^o) can conveniently be calculated numerically by the Lanczos continued 
fraction method 23 or kernel polynomial expansion 24 . We found essentially identical results for both 
algorithms. 16 However, the implementation of the Lanczos method turned out to be slightly faster. 
Well converged results were obtained already after several hundred iterations. 

We checked numerically that in the limit of large V, it is sufficient to calculate the spectral 
functions within the minimal Hilbert space. This enables us to study rather large systems and to 
address the question whether or not the two defects created by injecting one particle are closely 
bound to one another or not. 

The concept of a spatial separation of fcp's can be confirmed by comparing results within the 
minimal Hilbert space with those from an artificially restricted calculation in which a particle is 
prevented from decaying into two defects with charge e/2 each. In our finite-cluster calculations, 
a broad low-energy continuum is seen in the spectral functions for the unrestricted case, which is 
missing when the restriction is imposed. 22 The respective bandwidths are about 13t and St. This 
suggests a simple interpretation: The dynamics of two separated fcp's having a bandwidth of ~ 6t 
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(a) Spectral function A(k = tt/2,uj) for V = 25£ calculated for a y/32 x a/32 cluster for the effective t—g 
Hamiltonian with (a) g — 0.01 and (b) g — 1. A Lorentzian broadening 77 = O.lt is used. 



each due to six nearest neighbors would explain the calculated 13£, while a confined added particle 
has a much smaller bandwidth 8t, which is close to 6t. 

Further insight can be obtained if the amplitude of the ring exchange g=12t 3 /V 2 in the effective 
Hamiltonian is considered an independent parameter which is no longer restricted to the regime 
g<Ct as enforced by tcF. Figure PVT hi compares results for the "physical" regime corresponding 
to the previously considered parameters g ~ t 3 /V 2 = O.Olt with those for g = t. In the latter case, 
the broad continuum at the bottom of the spectral density vanishes and a sharp (5-peak evolves 
instead. The latter should be viewed as a Landau quasi-particle peak. This suggests the following 
interpretation: The ring exchange term leads for the undoped case to charge order, which in finite- 
cluster calculations is destroyed when a particle is added by the separation of the two fcp's. These 
are weakly bound to each other. If g<C£, the distance of the two fcp's is larger than the system 
size considered and thus the excitations seem to be deconfined. An artificially increased g leads to 
much stronger confinement and the diameter of the bounded pair becomes smaller or comparable 
to the system size. The pair acts at low energies as one entity. The electron as a (however strongly 
modified) "particle" is recovered and a finite weight of the quasi-particle peak results. These findings 
suggest that for the physical regime V/t « 10 quasi-particles with a spatial extent over more than 
hundred lattice sites are formed. The huge spatial extent is expected to lead to interesting effect. 
E.g., one may find with increasing doping concentration a transition from a "confined" phase to 
an "fcp plasma" phase with yet unknown properties when the average distance of bound fcp-fcp 
pairs falls below the diameter of a single bound fcp pair. Other open questions for future research 
include a thorough investigation of the 3D pyrochlore lattice. Even though the checkerboard lattice 
and pyrochlore lattice show many similarities, there exist certain differences between them, e.g., 
due to the higher spatial dimensionality ASi Here, differences of the two lattices arise in the U(l) 
gauge theory which describes the low energy excitations of the considered systems. The related 
compact electrodynamic is always confining in 2+1 dimensions while it allows for the existence 
of a deconfined phase in 3+1 dimension^ Systematic exact diagonalization studies as well as the 
application of Monte Carlo techniques will provide a deeper insight and we hope to confirm the 
expected existence of a deconfined phase in a 3D pyrochlore lattice. In addition, a natural extension 
of the model considered in this paper is the inclusion of spin. This leads to a more realistic model 
and could provide a better link to experiments. 
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